High  speed  imaging  chopper  polarimetry 

Roy  M.  Matchko*  and  Grant  R.  Gerhart” 

*Battelle  Scientific  Services  Program,  1307  W.Remuda  Way,  Payson,  AZ  85541,  928-472-3179, 

rmatchko@msn.com 

**US  Army  Tank- Automotive  Research  Development  and  Engineering  Center,  Warren,  MI  48397, 

585-574-8634,  grant.gerhart@us.army.mil 


1 


Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
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Abstract 


The  Stokes  method  provides  a  seminal  method  for  the  determination  of  a  state  of  polarization  of  a  beam 
of  light  using  measurable  quantities.  This  paper  presents  a  methodology  for  the  rapid  acquisition  of  the 
four  Stokes  parameters.  We  have  been  able  to  acquire  images  at  the  rate  of  500  images  per  second  using 
a  high  speed  video  camera.  Using  the  high  speed  video  camera,  an  optical  chopper,  a  computer  system 
containing  a  frame  grabber  and  software  to  coordinate  all  elements  of  the  system,  the  four  Stokes 
parameters  have  been  acquired  in  6  milliseconds.  Included  in  this  paper  is  a  description  of  all  elements 
of  the  system,  the  calibration  of  the  camera  and  polarization  elements,  the  design  of  the  optical  chopper 
polarimeter,  and  corrections  related  to  rotation  of  polarization  elements  during  exposure.  A  discussion 
of  errors  includes  variations  in  pixel  sensitivity,  ill-conditioning  of  the  Stokes  parameters,  aperture 
clipping,  quasi-monochromatic  approximation  and  variations  in  the  Stokes  parameters  during  data 
acquisition.  Validation  and  testing  of  the  methodology  includes  a  comparison  of  empirical  and 
theoretical  skylight  polarization  parameters  and  values  related  to  known,  rapidly  changing  polarization 
states  using  several  different  bandwidth  filters. 

Keywords:  polarization,  polarimetry,  Stokes  parameters,  image  acquisition,  image  processing,  skylight 
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Introduction 


The  Stokes  methodology  is  a  seminal  method  to  determine  the  polarization  state  of  light\  This 
proeedure  requires  four  independent  intensity  measurements  to  eharaeterize  the  most  general  ease  for 
elliptically  polarized  eleetromagnetie  waves.  Eaeh  measurement  eorresponds  to  one  of  four  distinct 
optical  filter  arrangements.  The  Stokes  parameters  in  this  paper  are  defined  as  So,  Si,  S2  and  S3. 

Temporal  and  spatial  registration  errors  are  two  major  sources  of  error  prevalent  in  all  4-Stokes 
parameter  measurement  systems.  They  may  occur  because  the  scene  changes  temporally  during  the  data 
acquisition  time  interval.  These  non-stationary  effects  can  alter  the  four  intensity  measurements  and 
cause  significant  errors  in  the  final  determination  of  the  Stokes  parameters.  These  effects  can  be 
important  in  scenes  where  changing  sun  position  or  cloud  conditions  produce  rapid  changes  in  ambient 
illumination.  In  laboratory  experiments,  variations  in  temperature,  pressure,  density  or  concentration 
variations  can  cause  similar  effects.  Spatial  registration  problems  occur  when  the  scene  shifts  in  angular 
orientation  relative  to  the  camera  during  the  image  recording  process.  For  example,  if  four  spatially 
adjacent  camera  systems  simultaneously  record  imagery  associated  with  a  different  polarization  filter, 
each  image  will  record  a  slightly  different  perspective  of  the  original  scene.  This  technique  would 
reduce  temporal  registration  errors,  but  would  potentially  increase  errors  due  to  spatial  misregistration. 
These  four  images  would  have  to  be  remapped  or  warped  in  order  to  correct  for  parallax  effects.  In 
single  camera  systems  background  vibration,  noise  and  unstable  platforms  can  also  produce  spatial 
misregistration  of  imagery.  Our  experience  has  been  that  these  effects  can  be  eliminated  or  controlled 
using  careful  experimental  procedures.  Temporal  registration  problems,  however,  are  typically  much 
more  difficult  to  correct. 
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Target  discrimination  and  object  recognition  are  major  objectives  in  our  pursuit  of  faster  imaging 
polarimeters.  Pseudo  ellipticity,  which  can  be  produced  by  temporal  registration  errors,  produces  false 
signatures  which  defeat  the  main  objectives.  The  authors  have  developed  an  imaging  chopper 
polarimeter  to  decrease  the  time  intervals  between  images  and  automate  the  image  acquisition  process. 
This  system  acquires  4-Stokes  parameter  imagery  using  a  high-speed,  digital,  video  camera.  It  reduces 
the  total  data  acquisition  time  from  30  seconds  in  our  earlier  COTS  digital  still  camera  polarimeter^  and 
from  0.6  seconds  using  our  rotating  retarder  apparatus  and  COTS  digital  camcorder  to  0.006  seconds  in 
our  new  imaging  chopper  polarimeter. 

The  Optical  Chopper  Polarimeter 

Figure  1  shows  the  Stokes  filter  wheel  used  to  acquire  the  four  Stokes  parameters.  It  consists  of  five 
trigger  holes;  three  linear  polarizers,  identified  by  1(0,0, 0),  1(0,90,0),  and  1(0,45,0);  one  circular 
polarizer,  identified  by  1(90,45,8)  and  an  opaque  reference  screen.  Figure  2  shows  the  components  of 
the  optical  chopper  polarimeter  and  the  connections  between  them. 

The  frequency  of  the  Stokes  filter  wheel  is  selected  from  the  optical  chopper  controller  (Scitec  300  CD) 
and  sent  to  the  Stokes  filter  wheel  via  connection  1.  The  computer  is  used  to  select  the  mode  of 
operation  (e.g.,  intemah  external  triggering,  etc.),  the  region  of  interest  (ROI,  the  pixel  area  to  be  used  in 
the  exposure),  and  the  exposure  time  of  the  video  camera  (Dalsa  1M75).  These  selections  are 
transmitted  to  the  video  camera  via  connection  4. 
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Signals  are  sent  from  the  infrared  sensor  to  the  optical  chopper  controller  via  connection  2  whenever  a 
trigger  hole  passes  between  the  infrared  transmitter  and  receiver.  The  optical  chopper  controller  sends  a 
signal  (e.g.,  a  5-volt  TTL  voltage)  to  the  frame  grabber  (BitFlow  R64)  via  connection  3  whenever  the 
optical  chopper  controller  receives  a  signal  from  the  infrared  sensor. 

The  frame  grabber  sends  a  triggering  command  to  the  video  camera  via  connection  4  whenever  the 
optical  chopper  receives  a  signal  from  the  infrared  sensor.  The  video  camera  sends  the  captured  image 
to  the  frame  grabber  via  connection  5  when  the  shutter  of  the  video  camera  closes.  Computer  software 
(StreamPix)  formats  (e.g.,  bmp,  jpg,  tif,  etc.)  the  raw  data  from  the  frame  grabber  and  saves  it  to  ram 
memory  or  to  the  hard  drive  of  the  computer. 

Calibration 

Calibration  of  the  Video  Camera 

Because  the  Stokes  parameters  use  intensity  measurements  or  a  parameter  proportional  to  intensity  and 
the  output  from  the  video  camera  is  in  RGB  gray  scale  values,  a  relationship  must  be  obtained  between 
the  RGB  gray  scale  values  and  the  intensity.  A  calibration  method  for  obtaining  this  relationship  is  to 
pass  a  beam  of  unpolarized  light  through  neutral  density  filters  of  different  optical  densities.  Curve¬ 
fitting  yields  RGB  gray  scale  values  as  a  function  of  optical  density  (OD).  The  best  curve-fit  was 
obtained  using  the  general  quadratic  equation 

A  +B  xy  +  C  y^  +  D  X  +  E  y  +  F  =  0  (1) 
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Using  X  to  represent  OD  and  y  to  represent  RGB  gray  scale  values,  Eq.  (2)  gives  OD  as  a  function  of 
RGB  gray  scale  values,  where  A,  B,  C,  D,  E  and  E  are  the  coefficients  of  Eq.  (1). 


-  (By  +  D)  +  7(By  +  D)"  -  4A  (Cy"  +  Ey  +  F) 
2A 


(2) 


The  calibration  curve  obtained  for  the  Dalsa  1M75,  when  unpolarized  incident  light  was  transmitted 
through  a  green  Kodak  #58  broadband  (100  nm)  filter,  was  a  hyperbola.  Relative  intensities,  I,  are 
obtained  from  optical  densities  as  given  in 


1=  10-OD 


(3) 


Calibration  of  the  Circular  Polarizer 

A  circular  polarizer  is  used  to  obtain  the  fourth  Stokes  parameter  and  is  identified  in  Fig.  2  as  1(90,45,8); 
8  is  the  phase  difference  of  the  two  orthogonal  components  Ex  and  Ey  of  the  electric  vector  transmitted 
through  the  circular  polarizer. 

A  circular  polarizer  consists  of  a  linear  polarizer  and  a  quarter-wave  retarder.  The  transmission  axis  of 
the  linear  polarizer  is  oriented  at  a  45°  angle  to  the  fast  axis  of  the  retarder.  The  phase  difference,  8,  can 
be  determined  empirically  by  positioning  a  linear  polarizer  in  front  of  the  circular  polarizer,  as  shown  in 
Fig.  3.  In  order  for  the  circular  polarizer  to  be  used  in  determining  the  fourth  Stokes  parameter,  the 
incident  light  must  first  be  transmitted  through  the  quarter-wave  retarder  and  then  through  the  linear 
polarizer.  The  empirical  8-values  are  obtained  when  a  =  0°  and  45°,  P  =  0°  and  0  =  -45°. 


6 


Using  Mueller  matriees"^,  the  intensity  of  the  light  transmitted  through  PI,  R  and  P2  of  Fig.  3, 1(a,P,0,8), 


is  obtained  from 


I(a,p,0,8)  =  lo  {cos^  (a+0-2P)  +  sin  2(a-P)  sin  (0-P)  cos  (0-P)  (1  +  cos  8  )}  (4) 


where  lo  is  the  intensity  of  the  incident  light  on  PI. 

Inserting  a  =  0°,  P  =  0°  and  0  =  -45°  and  a  =  45°,  P  =  0°  and  0  =  -45°  into  Eq.  (4)  and  solving  for  8 
yields 


s  =  acos  1  - 


1(45,0,-45,  g) 
1(0,0, -45, g) 


(5) 


To  obtain  empirical  8-values,  collimated  white  light  was  transmitted  through  a  frosted  glass  plate  and 
then  through  different  narrowband  (10  nm)  filters.  Using  a  Minolta  LS  100  luminance  meter  to  record 
data  and  Eq.  (5),  the  phase  difference  of  the  3M  circular  polarizer  as  a  function  of  wavelength.  A,,  is 


8  =  -2.52  X  10'°^  +  5.27  -  0.53 


(6) 


Eigure  4  shows  8  as  a  function  of  wavelength.  A,. 

Calibration  of  the  Stokes  Filter  Wheel 

The  four  polarization  elements  used  to  obtain  the  four  Stokes  parameters  must  have  the  same  optical 
density  or  calculations  of  these  parameters  must  incorporate  compensation  for  differences.  Using  the 
Dalsa  1M75  camera,  uniform,  unpolarized,  collimated  light  was  transmitted  through  each  element  of  the 
Stokes  filter  wheel  (Eig.  1)  and  the  optical  density  for  each  of  the  four  elements  was  obtained  from  Eq. 
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(2)  and  the  relative  intensities  from  Eq.  (3).  Eaeh  of  the  four  relative  intensities  is  divided  by  the  largest 
of  the  four  relative  intensities  to  obtain  a  compensation  ratio  (CR)  for  each  element.  Each  of  the 
empirical  relative  intensities  is  multiplied  by  its  associated  CR. 

Corrections  Due  to  Rotation  During  Exposure 

The  intensity  of  polarized  light  transmitted  first  through  a  quarter-wave  retarder  and  then  through  a 
linear  polarizer,  as  shown  in  Eig.  5,  can  be  obtained  through  the  use  of  Mueller  matrices"^  and  is  given  in 
Eq.  (7). 

I(Q,  0,  8)  =  0.5  {So  +  [cos  20  -  sin  2Q  sin  2(Q  -  0)  (1  -  cos  8)]  Si  + 

[sin  2Q  cos  2(Q  -  0)  -  cos  2Q  sin  2(Q  -  0)  cos  8  ]  S2  -  sin  2(Q  -  0)  sin  8  S3}  (7) 

I(Q,  0,  8)  is  the  relative  intensity  for  a  particular  set  of  angles  fl,  0,  and  8,  where  8  is  the  phase 
difference  produced  by  the  retarder;  So,  Si,  S2,  and  S3  are  the  Stokes  parameters  of  the  incident  polarized 
light. 

Three  of  the  apertures  in  the  Stokes  filter  wheel  are  fitted  with  linear  polarizers.  Since  retarders  are  not 
present  in  these  apertures,  fl  and  8  are  zero.  Equation  (8)  gives  the  general  equation  for  the  intensity  of 
the  light  after  transmission  through  any  of  the  three  linear  polarizers  in  the  Stokes  filter  wheel. 

Ilp  =  1(0,  0,  0)  =  0.5  [So  +  (cos  20)  Si  +  (sin  20)  S2]  (8) 
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The  fourth  chopper  blade  aperture  contains  a  left  circular  polarizer  with  its  quarter-wave  retarder  facing 
the  incident  light  and  its  linear  polarizer  adjacent  to  the  camera.  Since  the  fast  axis  of  the  retarder  is 
positioned  45  degrees  from  the  transmission  axis  of  the  linear  polarizer,  the  general  equation  for  the 
intensity,  I4,  of  the  light  transmitted  through  the  circular  polarizer  is  given  by 

I4  =  1(0  +  7t:/4,  0,  e)  =  0.5  [So  +  (cos  20  cos  s)  Si  +  (sin  20  cos  s)  S2  -  sin  8  S3]  (9) 

The  set  of  angles  associated  with  each  aperture  of  a  stationary  Stokes  filter  wheel,  when  associated 
trigger  holes  are  positioned  at  the  sensor,  is  given  in  Eq.  (10). 


(Q,  0,  8)1  =  (0,  0,  0),  (Q,  0,  8)2  =  (0,  90°,  0),  (Q,  0,  8)3  =  (0,  45°,  0), 

and  (Q,  0,  8)4  =  (90°,  45°,  8).  (10) 

Using  Eq.  (10)  in  Eqs.  (8)  and  (9)  and  solving  for  the  Stokes  parameters  yields 

So  =  Ii  +  l2  Si=Ii-l2  S2  =  2l3-So  S3  =  [So  +  cos  8  S2  -  2  Uj/sin  8  (11) 

However,  since  the  apertures  rotate  through  an  angle  5  during  the  exposure,  as  shown  in  Eig.  6,  the 
actual  intensities  recorded  by  the  camera  are  the  average  intensities,  represented  by  which  can  be 
determined  from  the  average  value  of  a  function. 

If  a  linear  polarizer  rotates  through  an  angle  5  during  the  acquisition  of  Ilp,  then  using  Eq.8,  the 
measured  value  of  Up,  corresponding  to  fi,  I2  or  I3,  is  obtained  from  Eq.  (12). 
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(12) 


I* 


LP 


o  6-\-5  o  9^d 

+  —  f  cos  26^  d6' +  ^-  fsin26'd6' 

5  \  ^  I 


Integration  of  Eq.  (12)  yields 

Plp  =  1/2  {So  +  (1/25)  [sin  2(0+5)  -  sin  2(0)]  Si  -  (1/25)  [eos  2(0+5)  -  eos  2(0)]  S2}  (13) 


If  the  cireular  polarizer  rotates  through  an  angle  5  during  the  aequisition  of  I4,  then  using  Eq.  9,  the 
measured  value  of  I4  is  obtained  from  Eq.  (14). 


I*  -1 
2 


So  + 


Sj  cos^ 
S 


(9+(5 

j  eos  29d9  + 


Sj  cos^ 
S 


0+S 


(14) 


Integration  of  Eq.  (14)  yields 


1*4=  V2  {So  +  [cos  s/25]  [sin  (20  +  25)  -  sin  (20)]  Si  -  [cos  s/25]  [cos  (20  +  25)  -  cos  (20)]  S2  -  sin  s  S3 } 


Using  Eq.  (10)  in  Eqs.  (13)  and  (15)  yields 


(15) 


1*1  =  I*(0,0,0,5)  =  V2  [So  +  (sin  5/5)  (eos  5  Si  +  sin  5  S2)] 

1*2  =  I*(0,90,0,5)  =  1/2  [So  -  (sin  5/5)  (eos  5  Si  +  sin  5  S2)] 

1*3  =  I*(0,45,0,5)  =  V2  [So  -  (sin  5/5)  (sin  5  Si  -  eos  5  S2)] 

I4*  =  I*(90,45,s,5)  =  V2  [So  -  (cos  s  sin  5/5)  (sin  5  Si  -  cos  5  S2)  -  sin  8  S3]  (16) 

Solving  Eq.  (16)  for  the  four  Stokes  parameters  yields 
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So  =  (1*1  +1*2) 


51  =  (5/sin  5)  [(sin  5  +  cos  5)  I*i  +  (sin  5  -  cos  5)  1*2  -  2  sin  5  1*3] 

52  =  (5/sin  5)  [(sin  5  -  cos  5)  I*i  -  (sin  5  +  cos  5)  1*2  +  2  cos  5  1*3] 

53  =  [So  +  cos  8  (2  1*3  -  So)  -  2  l4*]/sin  8  (17) 


Note  that  both  So  and  S3  are  independent  of  5,  henee,  So  and  S3  are  rotationally  invariant!  Also,  as  5 
approaehes  zero,  5/sin  5  approaehes  one  and  Eq.  (17)  reduees  to  Eq.  (1 1). 


The  degree  of  polarization,  P,  the  polarization  azimuth  angle,  \\f,  and  the  polarization  ellipticity  angle,  %  , 
given  in  Eq.  (18),  are  obtained  from  the  Stokes  parameters  given  in  Eq.  (17). 


^Si'  +  S^  +  Sj 


tan  2y/  =-^ 


sin  2j: 


7Si"  +  S/+S3 


(18) 


Aperture  Clipping 

Aperture  elipping  oeeurs  when  part  of  an  aperture  of  the  Stokes  filter  moves  out  of  the  field  of  view  of 
the  eamera  during  exposure.  This  phenomenon  is  especially  prominent  at  high  image  rates  when  the 
exposure  time  is  too  long. 

Eigure  7  shows  the  dimensions  of  an  aperture  of  the  Stokes  filter  wheel.  Aperture  elipping  begins  when 
the  angular  rotation  of  the  ehopper  blade,  5,  exeeeds  18°  during  exposure.  During  exposure,  the 
maximum  amount  of  rotation  that  ean  oceur  for  a  single  pixel  is  approximately  65°.  Sinee  the  elements 
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of  the  Stokes  filter  wheel  rotate  around  a  eirele  with  eonstant  speed,  and  5  images  are  aequired  in  one 
revolution  of  the  ehopper  blade 

|5|  =  2  71  ft  (rad)  =  72  Ft  (deg)  (19) 

where  F  is  the  eamera  image  rate  and  t  is  the  exposure  time  in  seeonds.  The  maximum  image  rate 
attainable  with  the  Seitee  300  CD,  operating  with  5  apertures,  is  500  images  per  seeond.  The  angular 
rotation  of  the  Stokes  fdter  wheel,  using  500  images  per  seeond  and  an  exposure  time  of  0.5  ms,  is  18°; 
just  avoiding  aperture  clipping. 

Errors 

Photosite  Sensitivity  Variations 

Although  flat  frames  help  in  the  correction  of  photosite  sensitivity  variations  in  the  Dalsa  1M75  video 
camera,  the  sensitivity  of  each  pixel  can  vary  with  intensity  by  as  much  as  10%,  as  shown  in  Fig.  10. 
Figure  10  gives  the  sensitivity  ratio  for  individual  pixels  in  selected  rows  and  column  #40  of  the  CMOS 
image  sensor  for  different  incident  light  intensities.  The  sensitivity  ratio  is  the  ratio  of  a  pixel  value  in 
the  flat  frame  and  the  average  pixel  value  of  the  entire  fiat  frame.  The  data  associated  with  Fig.  10  was 
obtained  by  transmitting  a  uniform,  collimated,  white  light  beam  through  different  neutral  density  filters. 
Fluctuations  in  pixel  sensitivity,  camera  noise,  line  voltage  fluctuations,  mechanical  vibrations  and 
image  rate  all  contribute  to  errors  in  the  calculation  of  the  Stokes  parameters. 
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Ill-conditioning  of  Stokes  Parameters 

The  Stokes  parameters  are  prone  to  giving  false  polarization  parameters  when  low  intensities  are 
involved  in  their  calculation  or  when  they  are  near  or  equal  to  zero;  the  errors  in  calculating  the  P,  \|/  or  y 
values  can  be  very  significant.  \|/-values  are  undefined  for  circularly  polarized  light  and  are  ill- 
conditioned  for  y-values  near  ±458. 

Conversely,  errors  in  the  y-values  can  be  very  small,  especially  for  high  degrees  of  polarization  (P 
1.00).  According  to  Eq.  (18),  when  P  =  1.00,  y-values  depend  solely  on  S3,  which  is  independent  of 
rotation  of  the  chopper  blade  during  exposure  and  hence  is  rotationally  invariant. 

Quasi-Monochromatic  Approximation 

Light  is  considered  to  be  quasi-monochromatic  when  the  wavelength  range  is  small  compared  to  the 
mean  wavelength.^  The  set  of  four  Stokes  parameters  given  in  Eqs.  (11)  and  (17)  are  generally 
considered  to  be  valid  for  only  quasi-monochromatic  light.^  Eigure  8  shows  the  relative  transmission 
through  three  different  filters  as  a  function  of  the  incident  wavelength.  The  primary  purpose  of  the  “hot 
mirror”  is  to  eliminate  the  transmission  of  infrared  radiation.  The  Kodak  #58  filter  is  a  broadband  green 
filter  with  a  bandwidth  of  approximately  100  nm.  The  Oriel  #56561  filter  is  a  narrowband  green  filter 
with  a  bandwidth  of  approximately  10  nm.  The  Oriel  #56561  filter  is  the  only  filter  shown  in  Eig.  8  that 
satisfies  the  quasi-monochromatic  condition.  However,  it  will  be  shown  in  the  studies  that  follow  that 
the  use  of  the  broadband  Kodak  #58  filter  and  the  “hot  mirror”  filter  produce  similar  polarization 
parameters  as  obtained  from  the  narrowband  Oriel  #56561  filter. 
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The  phase  differenee,  8,  of  the  eireular  polarizer  used  in  the  Stokes  filter  wheel  is  wavelength  dependent, 
as  Eq.  (6)  and  Fig.  4  illustrate;  hence,  the  fourth  Stokes  parameter,  expressed  in  Eqs.  (11)  and  (17),  is 
also  wavelength  dependent.  Since  Eq.  (6)  can  only  accommodate  one  specific  input  wavelength,  errors 
are  introduced  into  polarization  parameters  dependent  on  S3;  namely,  P  and  y.  The  polarization  azimuth 
angles  are  completely  independent  of  the  phase  difference  of  the  circular  polarizer. 

Although  8  and  A,  are  not  related  linearly  throughout  the  visible  spectrum,  they  are  nearly  so  within  a 
100  nm  bandwidth  (see  Fig.  4).  Therefore,  to  a  good  approximation,  the  wavelength  corresponding  to 
the  average  value  of  8  for  a  selected  100  nm  bandwidth  is  the  one  corresponding  to  the  midpoint  of  the 
bandwidth  . 

Three  separate  studies,  using  each  of  the  three  filters  shown  in  Fig.  8,  were  used  to  test  and  validate  the 
quasi-monochromatic  approximation.  Each  study  collected  500  images  from  the  Dalsa  1M75  video 
camera.  Since  the  image  rate  used  for  each  study  was  500  images/sec  and  five  consecutive  images  are 
required  to  determine  the  Stokes  parameters  (four  polarization  images  plus  a  blank  reference  image), 
100  different  polarization  forms  were  used  in  the  analysis  of  each  study. 

An  exposure  time  of  0.5  ms  was  used  to  avoid  aperture  clipping  during  exposure.  According  to  Eq.  (19) 
and  Fig.  7,  the  combination  of  500  images/sec  and  an  exposure  time  of  0.5  ms  produces  18  degrees  of 
rotation  for  each  filter  wheel  aperture  during  exposure;  hence  aperture  clipping  is  avoided.  Matlab 
scripts  are  used  to  calculate  the  empirical  polarization  parameters.  Equation  (17)  is  used  to  correct  for 
rotation  during  exposure. 
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A  region  of  interest  (ROI)  of  256  X  256  pixels  was  required  to  aehieve  500  images/see  from  the  Dalsa 
1M75  video  eamera.  Eaeh  and  every  polarization  parameter,  in  all  three  studies,  was  obtained  from  the 
average  value  of  a  eropped  100  X  100  pixel  area. 

Rapidly  changing,  known  polarization  forms  were  obtained  from  a  fixed  linear  polarizer  and  a  rotating 
quarter-wave  retarder,  as  shown  in  Fig.  9.  In  the  present  studies,  the  linear  polarizer  is  fixed  in  a 
horizontal  position  and  the  quarter-wave  retarder  is  rotated  at  25  RPMs  with  a  precision  motor.  This 
configuration  produces  the  input  source  to  the  imaging  chopper  polarimeter. 

Equations  for  the  Stokes  parameters  incident  on  the  Stokes  filter  wheel  are  obtained  from  the  use  of  Fig. 
9  and  Mueller  matrices  and  are  given  in  Eq.  (20)  for  the  case  a  =  0. 

So  =  1  Si  =  cos^  2p  +  sin^  2p  cos  a  S2  =  sin  2p  cos  2p  (1  -  cos  a)  S3  =  sin  2p  sin  a  (20) 

Since  the  retarder  is  transparent  and  rotating  when  data  acquisition  begins,  a  method  must  be  described 
to  identify  the  position  of  the  fast  axis  of  the  retarder  when  the  first  image  is  captured.  The  method 
adopted  in  these  studies  uses  Eq.  (20)  with  p  =  0  to  identify  this  position.  When  p  =  0,  Si  =  1,  S2  =  0; 
and  S3  =  0.  Since  a  Matlab  script  generates  all  of  the  empirical  Stokes  parameters,  it  is  easy  to  identify 
that  image  containing  P  =  0.  Since  the  retarder  rotates  at  25  RPM  (150°/sec)  and  the  image  rate  is  500 
images/sec,  there  are  0.3  degrees  of  rotation  of  the  fast  axis  of  the  retarder  between  consecutive  images. 
Theoretical  polarization  parameters  P,  \|/  and  %  are  calculated  from  Eq.  (18). 
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Results  of  the  Quasi-Monochromatic  Study 


Figure  11  compares  the  empirical  and  theoretical  polarization  parameters  P,  \j/  and  x  for  each  of  the 
three  filters.  Figure  12  compares  the  differences  (A)  between  the  empirical  and  theoretical  polarization 
parameters  P,  \|/  and  x  by  superimposing  the  results  from  the  three  color  filters. 

The  Polarization  Azimuth  Angle 

The  polarization  azimuth  angle  (\|/)  is  calculated  from  Si  and  S2,  as  shown  in  Eq.  (18).  The  largest  \|/- 
value  errors  in  this  study,  as  seen  in  Fig.  12b,  occur  near  P  =  ±45°.  Referring  to  Eq.  (20),  the  theoretical 
value  for  S2  is  zero  when  P  =  ±45°;  accounting  for  the  larger  errors  found  in  Eig.  12b  for  these 
orientations  of  p.  Eigure  12b  clearly  shows  that  deviations  from  the  quasi-monochromatic 
approximation  do  not  affect  the  \|/-values,  they  are  independent  of  the  incident  wavelength.  Except  for 
the  case  of  P  =  ±45°,  Eig.  12b  shows  that  the  errors  in  the  \i/-values  are  less  than  4°  for  all  visible 
wavelengths  .  Note:  The  polarization  azimuth  angle  is  undefined  for  circular  polarization  (a  =  0°,  P  =  ± 
45°  and  8  =  90°). 

The  Polarization  Ellipticity  Angle 

The  polarization  ellipticity  angle  (x)  is  calculated  from  Eq.  (18).  Eor  normalized  Stokes  parameters, 
which  are  used  exclusively  throughout  this  study,  the  denominator  of  %  is  the  degree  of  polarization,  P. 
Since  the  incident  light  on  the  Stokes  filter  wheel  is  completely  polarized,  the  theoretical  P-value  is  P  = 
1.0.  Therefore,  the  theoretical  x-value  depends  solely  on  S3  and  is  independent  of  the  other  Stokes 
parameters.  Eor  example,  P-value  errors  have  little  effect  on  the  empirical  y-value.  This  is  so  because 
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the  empirical  P-values  range  from  1.0  to  only  1.1;  hence,  the  denominator  of  x  is  not  much  different 
than  unity  and  the  empirical  x-values  depend  primarily  on  S3.  Since  S3  is  rotational  invariant  as  the 
Stokes  filter  wheel  rotates,  x-value  errors  related  to  rotation  are  eliminated.  The  only  errors  remaining 
relate  to  deviations  from  the  quasi-monochromatic  approximation  and  errors  introduced  by  the  CMOS 
image  sensor.  With  the  exception  of  only  a  few  of  the  300  data  points,  the  x-value  errors  are  all  less 
than  4°,  as  shown  in  Fig.  12c.  Calculations  show  that  when  the  empirical  RGB  gray  scale  values  are 
changed  by  ±4  gray  values,  simulating  variation  in  sensor  sensitivity,  the  x-value  can  change  by  3°. 
Therefore,  if  small  variations  of  RGB  gray  scale  values  can  account  for  x-value  discrepancies  of  3°,  the 
errors  due  to  deviations  from  the  quasi-monochromatic  approximation  must  be  very  small. 

The  Degree  of  Polarization 

The  degree  of  polarization  (P)  is  calculated  from  Eq.  (18).  Since  So  =  1  for  normalized  Stokes 
parameters,  P  is  just  the  numerator  in  Eq.  (18).  Eigure  12a  shows  that,  with  the  exception  of  two  points 
out  of  the  300  plotted  points,  the  largest  error  in  the  P-value  is  10%;  the  difference  in  the  empirical  and 
theoretical  P-values,  for  the  majority  of  points,  is  6%.  The  largest  errors  occur  for  the  narrowband  filter, 
which  is  not  surprising.  The  narrowband  filter  has  the  largest  optical  density  and  requires  greater 
incident  light  intensity  than  the  other  filters  to  produce  the  same  RGB  values.  The  largest  intensity 
available  was  used  in  the  data  acquisition  related  to  the  narrowband  filter.  On  the  other  hand,  neutral 
density  filters  were  used  in  conjunction  with  the  other  two  filters  in  order  to  prevent  over  exposure. 

Diurnal  Skylight  Study 

The  imaging  chopper  polarimeter  was  used  to  acquire  diurnal  polarization  parameters  of  a  small  area  of 
skylight.  Using  an  image  rate  of  500  frames/s  and  an  exposure  time  of  0.46  ms,  empirical  values  of  P,  \\i 
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and  X  were  aequire  from  10:00  to  16:00  at  intervals  of  15  minutes,  during  elear  skies,  at  longitude  111° 
west  and  latitude  34°  north.  Empirieal  values  for  skylight  were  obtained  from  averaging  over  30  X  30 
pixel  areas  centered  northwest  at  an  altitude  of  5°.  Figure  13  gives  the  results  for  this  study. 


The  Degree  of  Polarization 

The  theoretical  skylight  polarization  parameters  are  obtained  using  the  Rayleigh  scattering  model. 
Equation  21  gives  a  modem  version  of  Ford  Rayleigh’s  theoretical  description  for  skylight  polarization^. 


P  = 


sin^Q 
1  +  cos^Q 


(21) 


Q  is  the  scattering  angle  in  the  observation  plane  containing  the  incident  solar  ray,  scattering  center  and 
observer  .  Lord  Rayleigh  assumed  the  incident,  unpolarized  light  is  scattered  by  small  isotropic 
particles  (e.g.  air  molecules  much  smaller  than  the  wavelength  of  light).  The  scattering  angle  Q  can  be 
expressed  in  the  more  useful  earth-horizon  coordinate  system,  using  spherical  coordinates,  as  follows: 


cos  Q  =  cos  a  cos  0  cos  (A  -  (p)  +  sin  a  sin  0 


(22) 


where  a  is  the  altitude  of  the  sun,  A  is  azimuth  of  the  sun,  0  is  the  altitude  of  the  point  of  observation  on 
the  celestial  sphere  and  (p  is  the  azimuth  of  the  point  of  observation.  We  use  the  convention  that  south 
corresponds  to  0°  and  west  to  90°.  The  difference  between  theoretical  and  empirical  P-values  The 
largest  difference  between  the  empirical  and  theoretical  \j/-values,  as  seen  in  Fig.  13a,  is  less  than  6% 
between  10:00  and  14:30;  a  maximum  difference  of  17%  is  acquired  at  15:00. 
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The  Polarization  Azimuth  Angle 

In  a  previous  paper^  we  derive  an  equation  for  the  polarization  azimuth  angle  of  skylight  in  the  earth- 
horizon  eoordinate  system,  as  given  Eq.  (23). 


sin  a  cos  6* -sin  6*  cos  a  cos(®-A) 

cos  y/  -  — , 

±yl-(sina:  sin  6*  -i-cos6’  cos  a  cos(^-A))^ 


(23) 


where  a  is  the  altitude  of  the  sun,  A  is  azimuth  of  the  sun,  0  is  the  altitude  of  the  point  of  observation  on 
the  eelestial  sphere  and  (p  is  the  azimuth  of  the  point  of  observation.  The  largest  differenee  between  the 
empirieal  and  theoretieal  \|/-values,  as  seen  in  Fig.  13b,  is  4.5°. 


The  Polarization  Ellipticity  Angle 

Aeeording  to  Rayleigh  seattering  theory,  the  elliptieity  angle,  y,  of  skylight  is  zero.  The  largest 
differenee  between  the  empirieal  and  theoretieal  y-values,  as  seen  in  Fig.  13e,  is  4.3°. 

Equations  (21)  and  (23)  do  not  take  into  aeeount  interaetions  between  the  ineident  solar  radiation  and 
large  atmospherie  partieles  (eloud  droplets,  iee  erystals,  aerosols,  ete.)  and  effeets  due  to  multiple 
seattering  and  moleeular  anisotropy. 


Conclusion 
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Current  use  of  the  optieal  ehopper  polarimeter  has  produeed  outstanding  laboratory  results.  The 


differenee  between  empirieal  and  theoretical  \|/-values  is  less  than  4°,  except  for  the  case  of  P  =  ±45°, 
which  produced  a  10°  difference.  The  difference  between  empirical  and  theoretical  x-values  is  less  than 
4°.  The  largest  error  in  the  P-values  is  approximately  10%;  however,  for  the  majority  of  points,  the 
errors  are  less  than  6%.  The  use  of  the  quasi-monochromatic  approximation  with  broadband  filters 
shows  negligible  effects  on  the  x-values  and  only  small  effects  on  the  P-values. 


Future  studies  will  attempt  to  further  reduce  the  acquisition  time  of  the  four  Stokes  parameters. 
Additional  studies  will  use  the  optical  chopper  polarimeter  in  daylight  scenes  to  study  non-uniform, 
natural  and  man-made  targets. 
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Figure  Captions 


Fig,  1  The  Stokes  filter  wheeel. 

Fig.  2  The  eomponents  of  the  optieal  ehopper  polarimeter  and  the  eonneetions  between  them. 

Fig,  3  Apparatus  for  obtaining  the  phase  differenee  of  a  eircular  polarizer. 

Fig.  4  The  phase  differenee  of  a  3M  eireular  polarizer  as  a  funetion  of  wavelength. 

Fig.  5  Ineident  light  transmitted  through  a  retarder  and  a  linear  polarizer. 

Fig.  6  Rotation  of  the  Stokes  filter  wheel  during  exposure.  If  the  polarization  elements  rotate 
cloekwise,  then  5  is  negative. 

Fig.  7  The  dimensions  of  a  Stokes  filter  wheel  aperture. 

Fig.  8  Transmission  as  a  funetion  of  wavelength  for  a  “hot  mirror”  filter,  a  Kodak  #58  filter  and  an 
Oriel  #56561  filter. 

Fig.  9  Components  for  produeing  known  polarization  parameters;  the  input  source  to  the  chopper 
polarimeter. 

Fig.  10  Variations  in  the  sensitivity  ratio. 

Fig.  11  Comparison  of  empirical  and  theoretical  polarization  parameters,  (a,  b,c)  Degree  of 
Polarization,  (d,  e,  f)  Azimuth  Angle,  (g,  h,  i)  Ellipticity  Angle.  Input  source  is  a  25  RPM  retarder  and  a 
fixed  horizontal  linear  polarizer  (as  shown  in  Fig.  10).  Chopper  exposure  time  is  0.5  sec.  at  500 
frames/sec. 

Fig.  12.  Comparison  of  the  differences  (A)  between  the  empirical  and  theoretical  polarization 
parameters  P,  \\f  and  y  by  superimposing  the  results  from  the  three  color  filters,  (a)  degree  of 
polarization,  P  (b)  polarization  azimuth  angle,  \|/  and  (c)  polarization  ellipticity  angle,  y.  Input  source  is 
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a  25  RPM  retarder  and  a  fixed  horizontal  linear  polarizer  (as  shown  in  Fig.  10).  Chopper  exposure  time 
is  0.5  see.  at  500  frames/sec. 

Fig.  13.  Comparison  of  empirieal  and  theoretieal  skylight  polarization  parameters  during  elear  skies, 
(a)  degree  of  polarization,  (b)  polarization  azimuth  angle,  (c)  polarization  ellipticity  angle.  Chopper 
exposure  time  is  0.5  see.  at  500  images/sec. 


24 


Fig.l 


25 


Circular  Linear 

Polarizer  ^  Polarizer 


Fig.  3 


27 
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Fig,  6 
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Fig.  7 
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Fig.  8 
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Fig.  9 
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Fig.  11 
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